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I. INTRODUCTION 



Quantum spin chains are in fact three-dimensional magnetic compounds in which the 
magnetic interaction in one direction dominates over the two others. They are not mere 
academic curiosities since many real-world compounds do have this property. One such 
compound is KCUF3 , a satisfactory realization of the one-dimensional spin 1/2 antiferromag- 
netic Heisenberg model (lj. Early description of its crystallographic properties |2| confirms 
that it comes in two types, (a) and (d) HQ, both with similar spin structures. The (an- 
tiferromagnetic) coupling constant along the chain direction is roughly 100 times stronger 
than the (ferromagnetic) interchain coupling constant, a feature confirmed by electron-spin 
resonance measurements. This justifies amply a dominant one-dimensional behavior, and 
additional specific heat measurements are consistent with this interpretation. Also, an ideal 
strict one-dimensional magnetic chain will not exhibit long-range order at any finite temper- 
ature T [5|, but realistic quasi-one- dimensional systems have a characteristic temperature 
below which there is long-range order. For example, in KCUF3 long-range order is manifest 
below 38K for type (a) and 20K for type (d) Q. 

Though presumably simpler than higher-dimensional systems, quantum spin chains show 
strong and rich quantum behavior. If for example we consider antiferromagnetic systems, 
we would classically anticipate a Neel state traversed by spin-waves, and (linear) spin-wave 
theory is traditionally the usual framework for such a description. In this context, the 
excitation spectrum for the one- dimensional antiferromagnetic Heisenberg model is predicted 
to be f|: 

u cl (k) = 2 |sinfc| , (1.1) 



where k is the momentum transfer along the chain direction and here, the magnetic coupling 
constant is normalized to unity for convenience. But the true ground state is actually 
different from the Neel state 7|, and the lowest-lying excited states satisfy the so-called 



des-Cloizeaux- Pearson (dCP) dispersion relation 

^dcp (k) = 7T I sin k\ . (1-2) 
At first, these excitations were understood as spin-wave-like states with spin one and early 



experiments on CPC 



ons were understood as sp 
Qjjiand KCuF 3 [n[ Q 



were interpreted confirmation of this. 



But it was later shown |l2l,Uj| that the natural excitations of the model have spin 1/2 and 
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hence are fermions. These elementary excitations are called spinons and come always 
in pair. Furthermore, if the spin of the system is an integer, the spinons are bound and lead 
to well-defined spin- wave-like modes exhibiting a Haldane gap |l6j , a feature not seen in the 
(linear) spin-wave theory. If the spin is a half-integer, the spinons are unbound and there 
is no Haldane gap. Spin 1/2 compounds are even more interesting in that the spectrum 
is not simply given by a definite dispersion relation and indeed, analytic [l^l and finite 
chain ^| calculations showed that it is actually a continuum of excitations confined to first 
approximation, for a given k, between a lower bound uji (k) and an upper bound uj u (k) such 
that: 

uji (k) = uj dCP {k) ; uj u (k) = 2tt |sin k/2\ . (1.3) 

The ground state properties and more particularly the excitation spectrum of a mag- 
netic system are analyzed using inelastic neutron scattering, the cross-section of which is 
a function of the energy and momentum transfers u and k respectively 

QQ- All the 

above mentioned features regarding the Heisenberg model are confirmed by experiments. 
For example, inelastic neutron scattering on KCUF3 shows significant contribution to the 
scattering from regions consistent with the spin-wave continuum (jl.3j) and not with linear 



spin-wave theory, most particularly at high frequencies 



2l|. This implies the inadequacy of 



spin- wave theory |22j at such high energies in the quantum limit s = 1/2, this even if the 



24] 



two-magnon interactions are taken into account beyond the linear approximation j23l . 
But there is consistency at low-energies, and one must mention that spin-wave theory works 
fine for compounds with high spins (classical limit) like KFeSo with s = 5/2 where it gives 



accurate peak positions, line shapes and relative intensities |21j. One should also note that 



inelastic neutron scattering on KCuF 3 showed consistency of the field theory approach 



25]. 



valid only at long wavelengths, particularly the temperature dependence of the scattering 
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. This is important because it is an additional confirmation of the fermionic nature 
of the elementary excitations, the spinons. 

The inelastic neutron scattering is theoretically analyzed with the help of the two-point 
dynamic structure factor (DSF) S, for, in the Born approximation which is amply suffi- 
cient for the present purposes, the differential cross-section per scattering solid angle Q and 



outgoing neutron energy Ef is [19, 20]: 



d 2 o N kf 
dVl dEf ki 



9 -F (k)| 2 (*« - kA) S ij («, k) . (1.4) 



N is the number of scatterers, o"o = 0.2896b a unit for magnetic scattering, g the Lande 
factor, F (k) the magnetic form factor, kf the outgoing and ki the incoming neutron momenta 
respectively, i and j the cartesian coordinates. The DSF S is the Fourier transform of the 
two-point spin correlation function: 

i r+oo 

S ^ k) = 2^hN j „ dt S e iMTl - rm) -" t] (SI (0) SI (t)) , (1.5) 

where Yu m \ is the position of the spin Su m ) on the chain. The Hamiltonian corresponding 
to the antiferromagnetic isotropic Heisenberg model commutes with the total spin operator, 
which implies that S l i (u, k) is diagonal in % and j and S u (u, k) is the same for i = x,y and 
z. The averaging in (jl.5j) is generally done at finite temperature, but at low temperatures, 
only the ground state is retained. Also, for a quasi-one-dimensional system, the cross-section 
for neutron scattering with momentum transfer k depends only on the component k of k 
parallel to the chain. 

Historically, the first attempts to calculate the DSF for the antiferromagnetic Heisenberg 
model were made in the context of spin-wave theory which, for the chain lying in the in- 



direction, gives the transverse response as 

S^ T (u, k) = S y s l T (u, k) = |tan k/2\ 5 {to - u cX (k)) , (1.6) 

and the longitudinal component Sg WT (u, k) with a logarithmic singularity at uj c \ (k). Non- 
isotropy here is a consequence of the assumption of a Neel ground state with long-range 
order. Now the true quantum ground state does not have long-range order, which means 
isotropy must apply. The inadequacy of (jl.6|) triggered further efforts. First those of [28], 
where, guided by exact results of the one- dimensional XY model, numerical calculations on 
finite chains and known sum rules, it is constructed an ansatz for the DSF of the spin 1/2 
chain at zero temperature, the so-called Miiller ansatz: 

sjU( M .*)= e( "~"' (t j ) e i;- w ~" ) , (1.7) 

y/u 2 - Ujf (k) 

where G is the Heaviside step function. This form for the DSF has two main features, 
(i) A square-root singularity at the lower boundary of the spin-wave continuum, (ii) A 
built-in strict restriction to the spin-wave continuum itself, which is not physical at the 



high-frequency end if the Miiller ansatz is to represent the total DSF. 
ansatz predicts an abrupt high-energy cutoff while computer calculations 



ndeed, the Miiller 
28j suggest there is 



contribution of states outside the continuum, particularly from above, but with a scattering 
weight two orders of magnitude lower than that of neighboring states inside the continuum. 
Within the spinon picture, the small contributions outside the continuum are identified with 
processes in which more than two spinons are created. There is also a relatively small but 



systematic underestimation of the relative spectral weight near k = 7r |2jJ. But in overall, 
the Miiller ansatz gives reasonable results (to order unity) for known sum rules, is consistent 
with quantum Monte Carlo calculations |3] and supported by j^. More importantly, it 
is in good agreement with the inelastic neutron scattering experiments done on KCuF 3 
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The second set of efforts was to map the antiferromagnetic Heisenberg model in the 
long wavelength limit onto a relativistic quantum field theory and exploit the bosonization 
method of Luther and Peschel jsij to obtain an analytic expression of the DSF at finite 
temperature 25 . For integer spins, one finds that the DSF has a single mode with a gap 
as predicted in lfij |. For half-integer spins, the expression of the DSF at zero temperature 
agrees with the low-momentum limit of the Miiller ansatz and at finite temperature, there 
is good agreement with the measurements on KCUF3 j3, [^J . 

But all the above results regarding the dynamic structure function, though useful in their 
own right, are only approximate. This is important to note because over the decades, a 
number of quantum spin chains, most notably the Heisenberg model, have been amenable 
to exact solutions. As to the methods used, first there was the Bethe ansatz []| which, as 
already mentioned, gives the exact ground state for the model with characteristics differ- 
ent from the classical Neel antiferromagnetic ordering. Then techniques were developed in 
successful attempts to calculate exactly a number of thermodynamic quantities, techniques 
like the method of (commuting) transfer matrices and the Yang-Baxter equation j^|. These 
methods culminated in the so-called Quantum Inverse Scattering Method j^, a milestone 
towards the recognition of the quantum group symmetry present in the model. But up to 
here, all the exact results were for the most part concerned with only the static (thermody- 
namic) properties of the quantum spin chains; correlation functions that encode the dynam- 
ics remained elusive to exact treatment. One had to wait for advances in two-dimensional 
conformal field theory to see how the infinite two-dimensional conformal symmetry allows 
the computation of correlation functions of the so-called vertex operators, using bosoniza- 
tion techniques jsjj. The same strategy was then applied to quantum spin chains once it 
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became clear how to adapt the bosonization method to deformed commutation relations 



between creation and annihilation operators 



34, 



35 



36 



37| |. For the Heisenberg model, ac- 



tual formal manipulations are made in the context of the equivalent six-vertex model where 
vertex operators are defined in the framework of the infinite-dimensional representation of 
the quantum group. Correlation functions are then defined and, when mapped back to local 
spin operators, compact expressions of form factors are obtained [3| . 

The above treatment, though technically involved, opened the way for a systematic exact 
ireatment of the dynamic structure factor. In an expansion in the (even) number of spinons 
39I ] , it was first obtained an exact expression for the two-spinon contribution S2 j^J , and a 
com pari son with the Miiller ansatz showed that it gives a better account of the phenomenol- 
4^ [.Next an exact expression of the four-spinon contribution 5*4 to the DSF was 
43j and its behavior as a function of the neutron energy and momentum transfers 



ogy 
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derived in 

uj and k respectively was studied in 

The present work aims at furthering the investigation of the properties of the four-spinon 
DSF S4. Through an estimation of a number of sum rules, we address the issue of the 
spectral weight within S4 in the total S. The sum rules we use are: the static susceptibility, 
the integrated intensity, the total integrated intensity, the first frequency moment and the 
nearest-neighbor correlation function. We reach the conclusion that S4 contributes a weight 
between 1% and 2.5%, depending on the sum rule used, whereas S 2 contributes a weight 
between 70% and 75%. This is consistent with the finite chain calculations of j3] where 
it is observed that states just up the spin-wave continuum contribute with a scattering 
weight two orders of magnitude smaller than states within the continuum. The sum rules 
we investigate involve multi-dimensional integrations and these are performed using Monte 
Carlo algorithms. The statistical errors are reasonable. 

This article is organized as follows. After this introduction, we describe briefly in the 
next section the s = 1/2 antiferromagnetic Heisenberg quantum spin chain in the framework 
of give the definition of the dynamic structure factor 5* and decompose it into n-spinon 
contributions S n with n even. Then we write the expressions of S2 and S4 and give a brief 
account of their respective features and behaviors. In section three we describe the five sum 
rules we use to estimate the scattering weight of the four-spinon contribution and the results 
obtained. As already mentioned, the calculations use Monte Carlo integration methods and 
a discussion of the errors is given. Section four includes concluding remarks and indicates 
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few directions in which one can carry forward. 



II. HEISENBERG CHAIN AND DYNAMIC STRUCTURE FACTOR 

The antiferromagnetic s = 1/2 XXX Heisenberg chain is defined as the isotropic limit of 

the XXZ anisotropic Heisenberg model with the following Hamiltonian: 

+00 

H = "2 E A« + i) • (2-1) 

n=— 00 

A = (g + q~ 1 )/2 is the anisotropy parameter and the isotropic antiferromagnetic limit 
is obtained as A — > — 1~, or equivalently g — > — l~. Here cr^** are the usual Pauli 
matrices acting at the site n of the chain. Note that the definition (12. lj) and its 
isotropic limit are totally equivalent to the usual one if we consider the trans- 
formation U = exp (in J J^jS-oo s j)> which transforms the Hamiltonian H(J,A,h) = 
JEn^oc + «+i + A + <) into H(J,-A,h), i.e., f/tf ( J, A, />) f/- 1 = 

—H(J,—A,h). As mentioned already in the introductory section, we take the coupling 
constant J = 1 as well as H = 1 and the external magnetic field /i = 0. Lattice spacing is 
also taken equal to one. All this is consistent with the notation and conventions used in |3J| 
which we follow closely. 

The full exploitation of the quantum affine algebra U^s^) symmetry of the model re- 
quires an exact diagonalization of the Hamiltonian directly in the thermodynamic limit. 
Because of two different boundary conditions on the infinite chain, there are two equivalent 
vacuum states |0)j, i — 0, 1. The Hilbert space T consists of n-spinon energy eigenstates 

l£i,- ; Oei,-,en;* such that: 

n 

ff|^-,U« e »;i = E e fe)l^-'U^;i. ( 2 - 2 ) 

where e(£,-) is the energy of spinon j and £ ■ is a spectral parameter living on the unit circle. 
In the above relation, ej = ±1. The translation operator T which shifts the spin chain by 
one site acts on the energy eigenstates as follows: 

n 

T\£i, .-a., , = I! If 1> ""On .:i - , (2-3) 



i=l 



where r(£ 7 -) = e ip ^> and p(£ 7 - 



of e(£j) and are known [L3J, |43j and their isotropic limits are given in eq (J2-12) 
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, is the lattice momentum of spinon 7. The exact expressions 



below. The completeness relation in T reads: 



I = EE E h /ii^l: ^■■■>Ue 1 ,.,e n ;ii;e 1 ,...,eAl,---,a- i^A) 

j=i 



i=0,l n>0 {e ] =±l} j=lt „ 



A. The dynamic structure factor 



The dynamic structure factor we calculate is the zero-temperature limit of (|1.5j) up to 
27r, namely the Fourier transform of the transverse vacuum-to-vacuum two-point function 
defined by: 

/oo 
dtJ2e t{Mm) mai(t)a (0)10% , (2.5) 
■°° mGZ 

where u is the neutron energy, always positive, and k the neutron momentum component 
along the chain, denotes [a x ±ia y )/2. The DSF satisfies the following relations: 

S(uj, k) = S(uj, —k) = S(uj, k + 2tt) , (2.6) 

expressing reflection symmetry and periodicity. Inserting the completeness relation (J2.4|) 
and using the Heisenberg relation: 

a x m y ' z (t) = e iHt T~ m a*' v ' z (0) T m e~ iHt , (2.7) 

we can write the transverse DSF ()2.5|) as the sum of n-spinon contributions, with n even: 

S^-(u;,k)= S^-(uJ,k), (2.8) 

n even 

where the n-spinon DSF S n is given by: 



m£& ei,...,e n j 



x -v:;.:::,:^, a ^Uc-^i,-,-^) , (2.9) 

a relation in which X 1 denotes the form factor: 

K,.,JZl>-^n)= i<0k+ (0)|£ lf ...Ue,..,,^ (2.10) 

and z + m is to be read modulo 2. Note that each S n satisfies the symmetry relations (|2.6j) . 



The form factor X % is known 



38, 46]. It is expressed trace of vertex operators in 



the context of the infinite-dimensional representation of U^sL?). The trace is performed 



using g-deformed commutation relations of annihilation and creation operators the vertex 
operators are expressed with. From the form factor X 1 one can get a compact expression for 
the DSF S n in the general anisotropic case, an expression involving intricate complex 
contour inte gral s. The isotropic limit, the one of interest in this work, is obtained via the 
replacement : 

i = ie~ 2i£p , q = -e~ e , e -» + , (2.11) 

where p is the new spectral parameter suited for this case. The expressions of the energy e 
and momentum p in terms of p are: 

7T 

e (/°) = w 7 = — 71 s i n P '1 c °t P = sinh(27rp) ; — n < p < . (2.12) 

cosh(27rp) 

B. The two-spinon contribution 

The transverse two-spinon DSF S 2 is the less involved expression to derive from ([2.9)1 . It 
has been obtained in |41| and reads: 

- J W Q(u - u 2l (k))Q(u 2u (k) - u) 

The notation of the dynamic structure factor has been eased since we will deal only with 
the transverse DSF and the final results are independent of the vacuum state chosen. The 
function I(p) is given by: 

KHJ J t sinh(2t) cosh(t) ' v ; 

and (k) and uj 2 i {k) are the familiar upper and lower bounds of the spin-wave continuum 
of excitation energies given in ()1.3j) . The spectral parameter p is related to u and k by the 
relation: 

U 2u - U 2l 



S 2 (u, k - ix) = — — — ^ ^v^^ v^v (2 13) 



coshvrp=J^ f, (2.15) 

V U ~ U 2l 

which is obtained using eq ()2.12|) and the energy-momentum conservation laws: 

u = ei + e 2 ; k = -px - p 2 . (2.16) 



The properties of S 2 have been discussed in 



42 1 where a thorough comparison with the 



Miiller ansatz (|1.7|l is carried. We will simply note that: (i) The confinement of S2 in (|2.13|) 
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to the spin-wave continuum means that the Miiller ansatz (| 1.7(1 was in fact a description 
of the scattering weight from two-spinon processes only, and the spin-wave continuum is 
actually a 'two-spinon continuum', (ii) The exact two-spinon DSF S 2 exhibits a manifest 
square-root singularity at the upper boundary of the two-spinon continuum whereas the 
Miiller ansatz exhibits a square-root singularity at the lower boundary. 

C. The four-spinon contribution 

The expression of the four-spinon DSF 5*4 is given in 4^. For < k < tt it reads: 



S4 (u, k - 7r) = C 4 J dp 3 J dp 4 F (px, p 2 , p 3 , p A ) . (2.17) 



-7T — 7T 



For other values of k, it extends by symmetry using (|2.6(1 . The notation is as follows. C4 is 
a numerical constant equal to: 

2 -21 . -14 

C 4 = = 2.130627... x 10" 7 , (2.18) 

3|rQ)|Xf)| 8 

an expression in which T is Euler's gamma function and: 

f f 00 sinh 2 t (1 - 4-) expt\ 

A(z) = exp - / dt }, ™\ . 2.19 

K J \ Jo t sinh (2t) cosh t J v ; 

The integrand F in (|2.17(1 has a rather complicated expression. It writes: 

fto.ft.P^H E "-'-^■^^^'''"■ft-ft'^'. (2.20) 

(P1,P2) V « 

The different quantities involved in this expression are defined as follows: 

W = uj + 7i (sinp 3 + sinp 4 ) ; W u = 27r lsinX/21 ; K = k + p 3 + p 4 ; 

ccftpj = sinh (2^) , -Ti < pj < 0; (2.21) 
the function h in relation (|2.20j) reads: 

h( Pl ,p 2 ,p 3 ,p 4 )= J (%)> ( 2 - 22 ) 

l<«<i<4 
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where = p i — p- and the function / (p) is given in (|2.14|) . The function ge reads: 

4 

(2tt) 4 cosh faPj) 



./-J 

4 T 



x F — tt 1 TT 5 H (2-23) 

m _t^) n ^ sinh (ttp,,) T (m + 1 + 

In relation f)2.20|) . the sum is over the two pairs (pi,P2) and (p2,Pi) solutions of the following 
energy-momentum conservation laws: 

W = —7i (sinpi + sinp 2 ) ; K = —p\ — p2- (2.24) 

They read: 

( Pl ,p 2 ) = (-K/2 + arccos (W/ [2vr sin (K/2)]) , -K/2 - arccos (W/ [2vr sin {K/2)))) . 

(2.25) 

Note that the solution in (|2.25jl is allowed as long as Wi < W < W u where W u is given in 
(12~2TJ) and: 

Wi = TrlsinfTl . (2.26) 



The (analytic) behavior of the function F in fl2.2()j) was investigated in |43j ]. In particular, 
it was shown that the series gi is convergent for all values of its arguments and stays finite 
when two spectral parameters or more get equal. Since the function h goes to +oo in these 



regions j4lj, the integrand F of S4 is regular there. Furthermore, it was shown that F is 
exponentially convergent when one of the spectral parameters gets large, which means the 
double integration over and p^ in (12.17)1 is finite. All these analytic results pave the way 
for safe numerical manipulations. 

The behavior of S4 as a function of the energy and momentum transfers uj and k re- 
spectively is studied numerically in Q|. First is determined the extent in the (k, u;)-plane 
outside which 5*4 vanishes identically, i.e., the 'four-spinon continuum', by analogy with the 
two-spinon continuum. Then shapes of S4 as a function of u for different fixed values of k 
and vice versa are obtained. Consistency of S4 in three areas is obtained: (i) Confinement 
to the independently determined four-spinon continuum, (ii) Expected overall shape: sharp 
rise at the lower boundary of the continuum followed by a longer tail at the upper boundary, 
(iii) Similarity with the overall shape of the two-spinon contribution S%. 
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III. SUM RULES FOR THE DYNAMIC STRUCTURE FUNCTION 



ieisenberg antiferro- 
48 , 4t| . The ones we 



Once the expression of the four-spinon contribution to the total dynamic structure fac- 
tor is found and its behavior acquainted with, it is interesting to have an estimate of its 
contribution. One way is to consider sum rules, which are physical quantities related to the 
total dynamic structure factor, generally via specific integrals, the values of which we know 
exactly. This procedure is a reasonable good indicator of the weight of each contribution 
since each S n is positive, and so is their total sum S. 

For the present case, namely the one-dimensional nearest-neighbor 
magnet, a number of these sum rules has been derived in 47(, see also 
consider in this work are: the static susceptibility, the integrated intensity, the total inte- 
grated intensity, the first frequency moment and the nearest-neighbor correlation function. 
In the sequel, for each of these five sum rules, we calculate the corresponding contribution 
from the two-spinon DSF S2 and then the corresponding contribution from the four-spinon 
DSF S4. Then we make a weight comparison, with the corresponding exact result when 
possible, or simply between the two results. 

All the forthcoming calculations are numerical. The two-spinon calculations are fairly 
straightforward and standard quadratures are sufficient. But those related to S4 are quite 
more involved and we use a Monte Carlo algorithm for that. The x 2 factor and the standard 
deviation u 4 corresponding to each calculation are disp layed. Finally, note that all the exact 
results we refer to are known in the literature 



are disp l 

0,0 



A. Static susceptibility 

The first sum rule we consider is the one giving the static susceptibility k (k) for two local 

spin operators. It is related to the dynamic structure factor via the relation 19|: 

+00 

x(k) = - [ —S{u,k). (3.1) 
n J u 


The static susceptibility is known exactly in the limit k — > 0. It is found to be: 

x(0) = gmx(*) = ^a. (3.2) 

We therefore calculate the static susceptibility for both S% and S4 using (|3.1|) for small values 
of k, which we denote x 2 (k) and x 4 (k) respectively. Our results are displayed in Table 1. 
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Of course, x 2 (k) and X4 (k) for the chosen values of small k will not compare directly to 
x(0). Note also that it is not possible to take k directly equal to zero before integration 
because we would get zero identically. Finally, the behaviors of x 2 (k) and x 4 (k) themselves 
as functions of small k are not of interest in this study. 

One notes that the \ 2 factor stays around one but a 4/ x 4 is relatively high with respect to 
the coming sum rules, around 13% on average. This can easily be traced to the sensitivity of 
the integration to the low-energy region as there is the factor 1/u. Indeed, better precision 
is obtained when u appears in the numerator and not in the denominator. For this sum rule, 
X4 is only compared to x 2 and on average, the ratio X4/X2 is around 0.6%. We note that the 
estimated contribution of £4 with respect to this sum rule is smaller than the forthcoming 
ones. 



k 


x 2 (k) 


x 4 (k) 


04 


(T4/ x 4 


x 2 


x(0) 


x 4 /x 2 


0.1 


4.17340 10~ 2 


0.36927 10" 4 


0.6344 lO" 5 


17.1% 


0.813 


5.06606 10" 2 


0.1% 


0.2 


4.11826 10~ 2 


2.91744 10~ 4 


2.2114 10~ 5 


7.5% 


0.573 


5.06606 10" 2 


0.7% 


0.3 


4.00307 10" 2 


1.85373 10" 4 


1.8943 10- 5 


10.2% 


1.36 


5.06606 10" 2 


0.5% 


0.4 


4.48920 10~ 2 


3.99597 10" 4 


7.2076 10" 5 


18.0% 


0.897 


5.06606 10" 2 


0.9% 



Table 1: Static susceptibility for 62 and S4 



B. Integrated intensity 

The next sum rule we consider is the one that defines the integrated intensity with respect 
to the neutron energy transfer, namely: 

+00 

I (k) = — [ duS{u,k). (3.3) 

27T J 



The behavior of the integrated intensity is not exactly known for all values of k, but for 
small k, we know it is linear: 

I (jfc) ^ HA (k small) . (3.4) 

47T 

We will then use (13.3)1 to calculate, for small k, the integrated intensity I2 (k) coming the 
two-spinon DSF S 2 and the integrated intensity I4 (k) coming from the four-spinon DSF S4. 
The results are displayed in Table 2 and plotted in FIG 1. 
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We notice a clear linear behavior for h (k), with a slope roughly 75% of that of I (k). 
The four-spinon integrated intensity presents a somewhat linear behavior too, with a slope 
roughly 1% of that of / (k). The x 2 factor is still of order one and the standard deviation 
04 is here too about 13%, as for the static susceptibility. The explanation lies also in the 
small cu-region. 



k 


h(k) 


h(k) 


cr 4 




x 2 


I(k) 


hi I 


hi I 


0.1 


0.640435 10~ 2 


0.26462 10~ 4 


0.3656 10" 5 


13.8% 


1.16 


0.875535 10- 2 


73.1% 


0.3% 


0.2 


1.27924 10~ 2 


1.72237 10~ 4 


1.7739 10" 5 


10.2% 


1.30 


1.750704 10- 2 


80.2% 


1.1% 


0.3 


1.911885 10- 2 


2.98517 10~ 4 


4.9236 10~ 5 


16.3% 


0.607 


2.62606 lO" 2 


72.5% 


1.1% 


0.4 


2.643305 10~ 2 


3.11577 10~ 4 


4.0768 10" 5 


13.1% 


0.951 


3.501409 10" 2 


75.0% 


0.9% 



Table 2: Integrated intensity for S 2 and S4 



FIG. 1: Integrated intensity for S 2 and S4 for small k. 
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C. Total integrated intensity 



The total integrated intensity with respect to both energy and momentum is defined by: 

7T 7T OO 

I = - J dkl(k) = ^ J dk J duoS (to, k) . (3.5) 







It is known to be exactly equal to [4 



1 

4' 



(3.6) 



Table 3 displays our results for this sum rule. We note a better standard deviation: a 
relative value of 8% only, compared to the previous 13%. The \ 2 factor is of order one. The 
two-spinon contribution stays at around 74% whereas the four-spinon contribution is now 
at about 2.5%. This trend, namely a better precision and a 2.% contribution, is going to 
continue in the remaining sum rules. 



h 


h 


a 4 


a A /h 


x 2 


/ 


hi I 


hi I 


0.1845412 


0.00629916 


0.0004973 


7.9% 


0.430 


0.25 


73.8% 


2.5% 



Table 3: Total integrated intensity for S2 and S4 



D. First frequency moment 

The next sum rule we look at is related to the first frequency moment, defined by: 

00 

K(k) = — [ duouo S (u, k) . (3.7) 

27T J 







It is also known exactly, but for all k this time. It reads: 

K(k) = 8 -(ln2-\y^^j. (3.8) 

We have carried out the runs for small values of the momentum transfer and our results for 
the two-spinon contribution K 2 (k) and the four-spinon contribution K4 (k) are displayed in 
Table 4 below. They are also plotted in FIG 2. 

On sees that K2 (k) has the right quadratic behavior from ()3.8|) with a relative coefficient 
of roughly 69%. But what is remarkable is the fact that also K4 (k) has the same quadratic 
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behavior, with a stable coefficient of about 2.5%. The \ 2 factor is of order one and the 
standard deviation is better, around 6.7%. The better statistics comes partly from the fact 
that the neutron energy uj is present in the numerator and not in the denominator in the 
integral (|3.7j) . which makes the numerical integration more stable close to u — 0. One may 

oo 

then expect the higher frequency moments Ki(k) = J duo uo l S (u , k) with % > 1 to give 

n o 

even more stable results. 



k 


K 2 {k) 


K 4 (k) 


<7 4 


04/^4 


x 2 


K(k) 


K 2 /K 


K 4 /K 


0.1 


0.2053042 10~ 2 


0.80744 10~ 4 


0.6305 10" 5 


7.8% 


1.196 


0.2951853 10" 2 


69.5% 


2.7% 


0.2 


0.813574 10~ 2 


3.08253 10" 4 


2.1058 10" 5 


6.8% 


1.549 


1.177792 10~ 2 


69.1% 


2.6% 


0.3 


1.799619 lO- 2 


5.42950 10~ 4 


2.8636 10~ 5 


5.3% 


0.208 


2.639001 10~ 2 


68.2% 


2.1% 


0.4 


3.294491 10~ 2 


12.70855 10" 4 


8.8941 10" 5 


7.0% 


0.524 


4.664213 10~ 2 


70.6% 


2.7% 



Table 4: First frequency moment for S 2 and 64 



FIG. 2: First frequency moment for £2 and S4 for small values of k. 
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E. Nearest-neighbor correlation 



The last sum rule we consider is related to the nearest-neighbor correlation function 
defined by the following relations: 







/ z z \ Eg 



(3.9) 



Its value is also known exactly |47| : 



(3.10) 



where Eg is the ground-state energy per site which, in the thermodynamic limit, is given 
by: 

£o = -Mn2-i). (3.11) 

Our results for the nearest-neighbor correlation function are displayed in Table 5 below. The 
two-spinon contribution is about 72% and the four-spinon contribution is around 2.4%. The 
(relative) standard deviation is about 12% and the x 2 factor is of order one. 



( Cr n <T n+l) 2 


(«+i) 4 


<T 4 




x 2 




() 2 /(> 


() 4 /() 


-0.105594 


-0.0035309 


0.000421 


11.9% 


0.580 


-0.1477157 


71.5% 


2.4% 



Table 5: Nearest-neighbor correlation for S2 and £4 



IV. CONCLUSION AND OUTLOOK 

In this work, we have calculated five sum rules for the exact four-spinon dynamic structure 
factor 5*4 in the spin 1/2 Heisenberg antiferro magnetic quantum spin chain. The calculations 
were all numerical using Monte Carlo integration methods. The statistics is satisfactory: 
the x 2 factor is always around one and the standard deviation range between 7% and 
13% of the calculated quantity. It is better for sum rules in which the neutron momentum 
transfer u appears in the numerator since this make the integration more stable when close 
to zero. 

From this analysis, we conclude that the four-spinon DSF S4 contributes to the total 
dynamic structure factor S a factor between 1% and 2.5%, depending on the sum rule used, 

17 



whereas the two-spinon contribution is between 70% and 75%. These results are consistent 

28j | that the spectral weight of states 



with the expectation from finite chain calculations 
outside the two-spinon continuum is roughly two orders of magnitude smaller than the 
spectral weight of neighboring states within the continuum. 

There are five directions in which one may wish to carry forward. The first direction is 
to try to determine an expression for the six-spinon contribution Sq, study its behavior as 
a function of u and k and then determine its contribution to the total S. But one should 
know that, technically, matters may be more involved. 

The second direction to explore is the study of the dynamic structure factor in the 
anisotropic case. This is also of physical interest since perfect isotropy is only an ideal 
limit. The model is exactly solvable and we do have generic expressions for S n in the form 
of contour integrals in the spectral parameters' complex planes j^. The difficulty here is 
that the integrands involve much more complicated functions which are already present in 
S2, and one should expect intricate complexities in this more general case. 

The third direction to explore is to include a (small) external magnetic field. There are 
finit e-chain captions in this regard, H and more recently Hi! Bn« one has t o 
remember that the model in not exactly solvable in this case. One will then have to try 
small perturbations around the zero-field limit solution. 

The fourth direction is the finite-temperature case. Here too there are finite-chain 

n 

and field-theory results [25J. It is then certainly interesting to see the temperature effects 
on S2 and perhaps on S4. 

The fifth direction is to look into the situation of a spin-one chain. The model is still 
exactly solvable and, exploiting t he q uantum group symmetry, compact expressions for the 



form factors are available 15 



HaMano gap Q through these 1 oxact manipulations. 



53 



541 ] . One key issue in this regard is to try to recover the 
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